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Abstract 



We discuss black holes in an effective theory derived from a superstring 
model, which includes a dilaton field, a gauge field and the Gauss-Bonnet 
term. Assuming U(l) or SU(2) symmetry for the gauge field, we find 
four types of spherically symmetric solutions, i.e., a neutral, an electrically 
charged, a magnetically charged and a "colored" black hole, and discuss their 
thermodynamical properties and fate via the Hawking evaporation process. 
For neutral and electrically charged black holes, we find critical point and a 
singular end point. Below the mass corresponding to the critical point, no 
solution exists, while the curvature on the horizon diverges and a naked sin- 
gularity appears at the singular point. A cusp structure in the mass-entropy 
diagram is found at the critical point and black holes on the branch between 
the critical and singular points become unstable. For magnetically charged 
and "colored" black holes, the solution becomes singular just at the end point 
with a finite mass. Because the black hole temperature is always finite even at 
the critical point or the singular point, we may conclude that the evaporation 
process will not be stopped even at the critical point or the singular point, 
and the black hole will move to a dynamical evaporation phase or a naked 
singularity will appear. 
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I. INTRODUCTION 



One of the most fascinating dreams for all physicists is the unification of all fundamental 
forces, i.e., electromagnetic, weak, strong and gravitational interactions. The electromag- 
netic and weak interactions are successfully unified in the Weinberg-Salam theory. The 
strong interaction is described by quantum chromodynamics (QCD) and is likely to be uni- 
fied with the Weinberg-Salam theory into a grand unified theory in the context of gauge 
theory. The gravitational interaction, however, is not yet included, in spite of a great deal of 
effort. The most promising candidate for a unified theory of all interactions is a superstring 
theory, which may unify everything without any divergences. 

Although such a unified theory may become important in the strong gravity regime, 
however, little work has been done on extreme situations such as a black hole or the early 
universe. Since methods to study the strong gravity regime in a string theory are not well 
developed, most analysis have been performed by using an effective field theory inspired by a 
string theory, which contains the leading or next to leading order terms of the inverse string 
tension a'. One such application is string cosmology. Some puzzles in Einstein cosmology 
might be solved with a string theory. For example, while the singularity theorem demand 
that the universe has an initial singularity in the Einstein gravity, a string inspired model 
can remove it and provide a non-singular cosmology |]T],0]. 

Another application is black hole physics. The first study was made by Gibbons and 
one of the present authors in the Einstein-Maxwell-Dilaton (EMD) system || and the same 
solution was also discussed in a different coordinate system in ref. [|J. They found a static 
spherically symmetric black hole solution (GM-GHS solution) with dilaton hair. Since the 
dilaton hair cannot appear without electromagnetic hair, it is classified as secondary hair. 
After this work, many solutions were discussed in various models. Thus the following ques- 
tion naturally arises: How are the black hole solutions affected by the next to leading order 
term in a', in particular by the higher-curvature term? This was first considered inde- 
pendently of a string theory. Wheeler analyzed the effect of the Gauss-Bonnet (GB) term 
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without a dilaton field ||. When the dilaton field is absent, the GB term does not give 
any contribution in a four dimensional spacetime because it becomes a surface term and 
gives a topological invariant. Then he studied black holes in spacetime with more than four 
dimensions. These solutions are called dimensionally continuum black holes ||. Callan et 
al. |7|] discussed black hole solutions in the theory with a higher-curvature term R fJiUpa R IJ ' upa ' 
and a dilaton field, and Mignemi and Stewart |§ took both the GB term and a dilaton field 
into account in four dimensional spacetime. In their work, field variables are expanded by 
the inverse string tension a' and the first order terms of a' are taken into account. Using 
this perturbation, they constructed analytic solutions. There are also some studies which 
clarify the effect of an axion field as well as the dilaton field, a U(l) gauge field and the 



GB term P-pl- Either dyon solutions or axisymmetric stationary solutions are analyzed 
because the axion field becomes trivial if the gauge field does not have both electric and 
magnetic charges and the spacetime is static and spherically symmetric. In all these models 
containing a higher-curvature term, a perturbative approach was used. Assuming the GB 
higher curvature term, recently, Kanti et. al. calculated a neutral solution without such a 



perturbation and found interesting properties [O . More recently Alexeyev and Pomazanov 



also discussed the internal structure of these solutions JTB| ]. 

When we turn to another aspect of a unified theory, we find a non-Abelian gauge field. 
One of the most important facts about black holes with a non-Abelian gauge field is that 
we have the so-called colored black hole, which was found in Einstein- Yang-Mills system 
T4| soon after the discovery of a particle-like Bartnik-McKinnon(BM) solution |15|]. These 



solutions can exist by a balance between the attractive force of gravity and the repulsive 
force of the Yang-Mills (YM) field. If gravity is absent, such non-trivial structures cannot 
exist. In this sense, these solutions are of a new type. Although both the BM particle and 
the colored black hole are found to be unstable against radial linear perturbation |]IE|| , they 
showed us a new aspect of black hole physics and forced us to reconsider the black hole 
no-hair conjecture. 

After these solutions were discussed, a variety of self-gravitating structures and black 



hole solutions with a non-Abelian field were found in static spherically symmetric space- 
time p)-|T9|j. The Skyrmion [^0|,^lj or the Skyrme black hole in Einstein- Skyrme 



system, the particle solution with a massive Proca field or the Proca black hole [jzl] in 
Einstein-Proca system, the monopole p5|- p7| , p!9[1 or the black hole in a monopole [p6|- |28|JT9[] 
in the Einstein- Yang-Mills-Higgs (real triplet) system and the sphaleron [^,0 or sphaleron 
black hole [^4J] in Einstein- Yang-Mills-Higgs(complex doublet) system have been discovered. 
A particular class of the Skyrme black hole, the Proca black hole and the monopole black 
hole turns out to be stable against radial perturbations. In particular the monopole black 
hole can be a counterexample of the black hole no-hair conjecture |n|, because it is highly 



stable and is formed through the Hawking evaporation process from the Reissner-Nordstrom 
black hole. We also investigated the dilatonic BM particle and the dilatonic colored black 



hole solution in the Einstein- Yang-Mills-Dilaton (EYMD) system p3| , |30| , which are direct 
extensions of the GM-GHS solution. 

Although those non-Abelian fields may be expected in some unified theories, such black 
holes must be very small and then some other contributions such as the GB term and/or the 
moduli field may also play an important role in their structure, if the fundamental theory 
is described by a string model. Donets and Gal'tsov showed that a particle-like solution 
does not exist in the EYMD system with the GB term ||31|| . Then, they assume a numerical 
constant (3 in front of the GB term, where (3=1 corresponds to the effective sting theory. 
They showed that there is a critical value (5 cr = 0.37, beyond which no particle-like solution 
exists. However, we expect that a black hole solution can exist in this system, so we also 
study the case of a SU(2) YM field. 

In this paper, then, we study black holes in a theory inspired by a string theory, i.e., 
in a model with a dilaton field, a gauge field, and the GB curvature term, and discuss 
their properties. This paper is organized as follows. We outline the a model and field 
equations in section 2, and present various types of new solutions (neutral, electrically 
charged, magnetically charged and "colored" black holes) in section 3. In section 4 we study 
the thermodynamical properties of those black holes. Section 5 includes discussions and 



some remarks. 



II. A MODEL AND FIELD EQUATIONS 

We shall consider the model given by the action 



where k 2 = 8nG, and discuss a spherically symmetric, static solution. This type of action 
comes from low-energy limit of the heterotic string theory ||32|| . The dilaton field is (f) and 
7 = \/2 is the coupling constant of the dilaton field to the gauge field, if is a three form 
expressed as 

a' 

H = dB + — (n 3L - n 3Y ) , (2) 

OK 

where B^ v is the antisymmetric field in the gravitational multiplet. Q 3 l and fl 3 y are the 
Lorentz and gauge Chern-Simon terms respectively; 

fi 3L = tr (lu A R - A uj A cuj , (3) 
Q 3Y = Tr (a A F - A A A Aj , (4) 

where the traces are taken over the Lorentz and gauge indices. to afll/ = e^e^V p e I/ (b) is the 
spin connection with vierbein Using the dual of the Bianchi identity, we can rewrite a 
part of the action 

- le-^H^H^ = -^e 2 ^%b KR d% KR 



+ 



a' 



b KR eW (Ra^R^f + TrF pa F^)} . (5) 



The pseudoscalar field h^R is the Kalb-Ramond axion. If a spacetime is static and spherically 
symmetric and the gauge field does not have dyons but only an electric or a magnetic charge, 
which is the situation we will discuss here, then the second term of equation (|5]) vanishes. 
Then the axion field can be regarded as a massless scalar field. Such a scalar field, however, 



is trivial because of the black hole no hair theorem | 33| . Hence we can drop the axion field 
here. 

F is the field strength of the gauge field expressed by its potential A. If it is a U(l) 
gauge field (the electromagnetic field), then we consider only a electrically or magnetically 
charged black hole, i.e., 

A = adt and F = ^ dr A dt. (6) 

dr 

For SU(2) gauge field, we assume the Witten ansatz ||34|| , which is the most generic form of 
a spherically symmetric SU(2) YM potential. The YM potential becomes 

A = ar r dt + br r dr + [drg — (1 + w)t^\ d9 + [(1 + w)Tg + dr^\ sin 9d(p, (7) 

where a, b, d and w are functions of time and the radial coordinates, t and r. We have 
adopted the polar coordinate description (r r ,rg,r^), i.e. 

r r = -— [<r \ sin 6 cos + <r 2 sin 9 sin <fi + cr 3 cos 9], (8) 
Te = 2~-[ (T i cos ^ cos + °"2 cos 9 sin — cr 3 sin 9} , (9) 
= —[-tr! sin + cr 2 cos 0], (10) 

whose commutation relations are 

[r a ,r 6 ]=r c a, 6, c, = r, 9, or 0. (11) 

Here <Tj (i = 1, 2, 3) denote the Pauli spin matrices. We can eliminate b using a residual 
gauge freedom. In the static case, the part of the YM equations is integrated as d = Cw 
where C is an integration constant. We can set C = i.e. d = without loss of generality. 
The remaining functions a and w depend only on the radial coordinate r. As a result, we 
obtain a simplified spherically symmetric YM potential as 

A = a(r)r r dt - [1 + w(r)] T^dO + [1 + w(r)} r e sin 9d(f). (12) 

Substituting this into F = dA + A A A, we find the field strength: 



da da da 

F = —T r dr Adt + —Tsdr A dO + —Tgdr A sin 

dr dr dr 



1 — w ) T r d6 A sin Odd) + awTedt A dd + awr&dt A sin 



(13) 



Comparing (|T3|) with the field strength of the U(l) gauge field, we find that a and w play 
the roles of an electric and a magnetic potentials, respectively. This expression can be used 
for U(l) gauge field if we formally set r r = 1 and tq = t^ = 0. 
The GB term, R 2 , is defined by 



R 2 = R^R^ - AR^BT + R 2 . 



(14) 



This combination is introduced to cancel anomalies and has the advantage that the higher 
derivatives of metric functions do not appear in the field equations. Setting a'/n 2 = l/trg 2 , 
g is regarded as a gauge coupling constant. A numerical constant (3 is introduced in Ref. |3l| 
to find a non-trivial particle-like solution, but we fix it to be unity because we are interested 
in the effective string theory. 

Because of our ansatz, the metric is of the Schwarzschild type, 



ds 2 



2Gm\ 
I j e 



-»dt 2 + ( 1 - —) 1 dr 2 + r 2 (d6 2 + sin 2 9d<p 2 ). 



(15) 



The mass function m = m(r) and the lapse function 5 = 5(r) depend on only the radial 
coordinate r. 

Varying the action (|l|) and substituting ansatze (|12|) and fllB]) , we find the field equations; 
5' = h- 1 



1 e ~ 
-~rw 



4f 



^ ( <- 2m 

28 t}-2~2 2 , 12 i ■ LUl i 2 ±12 ±H\ 

e B a w + w H 7 <p — •yep ) 



(16) 



m! = h 1 



Bf 



-;2j/2 



-7</> 



4 

2m 



+ 



16 



(r~V 2 + 2B" 1 aW) + 2Bw' 2 + 
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2\2 



+45^ (7 V 2 - 70") + 87^ [B — - 



-5-2 d .1/ 

e r 



r 

/ p-T0 

7e r 



rn 



,25 J ~/2 



2dSv 2 



a + 



2«/ : 



2\2 ' 



-B 



(1-w 2 ) 

ji4 



+±{2f + B(S'f + f)} + ±(6'f-f) 



2e s e-^B~ 1 dw 2 = 0, 



e - s Be-^w'\ + e s e-^B~ l d' 2 w + e _6 e 
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s^ wjl- w 2 ) 



0. 



(17) 

(18) 
(19) 
(20) 



Here we have used the dimensionless variables; r = t/a/o 7 , m = Gm/\fa' and a = a\fa'. 
A prime in the field equations denotes the derivative with respect to f , and 

2m 



B = 1 



~ 5 

r 



-j4> 



m 



f = h-' 



ml ,9 e 



-70 



16r 



2<5 /~2~/2 OD -l~2 2\ oD '2 

e ( r a — IB aw- 2Bw 



[1-w 



2\2 



(21) 
(22) 
(23) 



Note that those equations can be applied for a U(l) gauge field by setting 



w = ±1, 



aw = 0, 



(24) 
(25) 



with a being a non-trivial potential. The latter condition fl25|) corresponds to the vanishing 
of the self-interaction due to the non-Abelian term. 

As for boundary conditions for the metric functions on the event horizon and at spatial 
infinity, we impose following three ansatze: 

(i) Asymptotic flatness at spatial infinity |35| , i.e., as r — ► oo, 

(26) 
(27) 



m(r) — > M = finite, 



6(r) -> 0. 



(ii) The existence of a regular horizon r#, i.e., 

2Gm H = r H , 
5h < oo. 

(iii) The nonexistence of singularities outside the event horizon, i.e., for r > r#, 

2Gm(r) < r. 



(28) 
(29) 



(30) 



The subscript H is used for the values at the event horizon r = r^. As for the field functions 
we have 



S 



0^0, (31) 

a -> 0, (32) 

!±1, (globally magnetically neutral solution). ^ 
0, (globally magnetically charged solution). 

as r — ► oo and impose their finiteness on the horizon. These conditions guarantee that 
the total energy of the present system is finite [[3(J. The boundary conditions of the field 



functions on the event horizon depend on the gauge group. Hence we will discuss them 
individually in the next section. 

III. NEW DILATONIC BLACK HOLES WITH GAUSS-BONNET TERM 

In this section we present the solutions of the field equations (|16|)~(|2"0"1). We classified 
them into four types by their gauge charge, i.e. neutral, electrically charged, magnetically 
charged and "colored" black holes. All solutions need numerical analysis. 

A. Neutral Black Hole 

First we consider the case without any gauge field, which is the simplest. The solution is 
a Schwarzschild type black hole modified by the GB term coupled to a dilaton field. Then 
the gauge field should be set to 

5 = 0, (34) 
w = ±l, (35) 



which satisfies the field equations (|T^) and (20) . From Eq. flT8|), we find the following 
relation for the dilaton field on the event horizon 

^ ~ % + 3 = °' (36) 
where A = e~ 1( ^ H / Af 2 H . The quadratic equation ( p6|) has two roots as 
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Hence, for each value of (j>H at a fixed event horizon f#, we have two possible boundary 



values <j)'±. 

We integrate the field equations (|16D~(18D from the horizon r = r# with the boundary 
conditions fl28|)~(|30D and Since the equation of the dilaton field (|18D becomes singular 
on the event horizon, we expand the equations and variables by power series of f — f # to 
guarantee the regularity at the horizon, and use their analytic solutions for the first step of 
integration. We show the behavior of the field functions of neutral black holes with three 
different radius of event horizon in Fig. 1. We find black hole solutions with a regular 
horizon only when we choose <p' H = <p' + . 

For smaller black holes, the dilaton field varies more rapidly than that of the larger 
ones. This means that stringy effects become more important for a smaller black hole, as we 
expected. We see that the mass function decreases first near the horizon and then increases 
afterward, approaching a finite value. There is a region where the effective energy density, 
which is defined by (dm/dr)/4:7ir 2 , becomes negative. From Eq. (|17D , rh' is evaluated on 
the event horizon as 

m » = " 2(1-^) = (38) 
which is negative definite. Hence the function m always decreases in the vicinity of the 

event horizon. It is one of the essential points for existence of the neutral black holes 

that the energy density becomes negative [0. Regarding the GB term as a source term 

of the Einstein equations and applying the similar analysis that was used to prove the 

no-hair theorem for a scalar field, we can show that if the time-time component of the 

energy-momentum tensor (the effective energy density) is positive everywhere, no non-trivial 

solution can exist. However in our situation this is not the case. So we find a new solution 

even if the gauge field is absent. 

We show M — th relations in Fig. 2. Note that there is an end point for each branch, 

where <p' + and <f>'_ coincide. With this fact, we can prove that <\)" H and S' H diverge. We have 
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also shown that R VLVpa W' vpa diverges on the event horizon, which means a naked singularity 



appears at the end point of the branch ||12|| . We shall call this point the singular point S. 

Near the singular point, we also find a critical point C, which gives a lower bound for the 
black hole mass, i.e., below which mass no solution exists (see Fig. 3(b)). Since this solution 
is a modification of Schwarzschild black hole, we find no other black hole solutions without a 
gauge field. No naked singularity appears for the solution at the critical point. We have two 
black hole solutions in the mass range of Mq < M < Ms, where Mc = 6.02771Mpi/g and 
Ms = 6.02813Mpi/g are the masses of the critical and singular solutions, respectively. As 
we will see later, the stability of these black holes changes at the critical point C indicates 
that the singular end point solution becomes unstable. When we discuss the evolution of 
the black holes, it will also important be that there exists a regular critical solution whose 
mass is smaller than that of the singular solution 



B. Electrically Charged Black Hole 

Electrically charged black hole solutions with the U(l) gauge field can be obtained by 
setting w = ±1, aw = in Eqs. (^6|) ^(|20| ) . Then the electric potential a" is integrated once 
to give 



,Qe 



~a> = e"V^, 



(39) 



where Q e is a constant of integration and denotes a normalized electric charge. The physical 
charge is given as gQ e . 

From Eq. flllf ), we find the following relation for the dilaton field on the event horizon 

<lA^ A e^»Ql$l - [A 2 je 2 ^ H Q 2 e (2A 2 7 2 + 5A 7 2 + l) - 2Ay 



+ 
+ 



A^e^Ql + Ae 2 ^«Q 2 e (2A 2 7 2 + 4A 7 2 + l) - 2 



^l e ^ H Ql - A^e 2 ^ H Ql (6 A + 1) + 6A 7 



0. 



(40) 



There are three roots of (f)' H , (/)[ < <p' 2 < 03- We can, however, obtain a regular solution 
only for (f)' H = <ft' 2 . This is understood from the fact that in the limit of Q e — > 0, we recover 
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the previous condition ( |36D and (p' 3 = oo. We plot the field functions of the solutions with 
Q e = 1.0 in Fig. 3. Their behaviors are almost same as those of the neutral case qualitatively. 

We show M — Th relations of black holes with Q e = 0.4 and 1.0 in Fig. 2. As the 
neutral black hole case, we also find the critical and the singular points (the end point on 
the branch), C and S. Fig. 2(c) is a magnification around the critical point C. No singular 
behavior appears for the critical solution and below the critical point, no solution exists. At 
the singular point, a naked singularity appears. 

The existence of a critical mass is also known in a Reissner-Nordstrom black hole with a 
fixed charge. In that case the outer and inner horizons coincide and the black hole becomes 
extreme at the critical mass. Our solution curve with a constant charge in M — r^ diagram 
is also vertical at the critical point. However, we will see in the next section that our 
critical point C has different thermodynamical properties from that of the extreme point in 
Reissner-Nordstrom black hole. 



C. Magnetrically Charged Black Hole 

Next we turn to a magnetically charged solution. In the present subsection we extend 
the gauge field from the U(l) to SU(2) group. This can be done by setting 

5 = 0, (41) 
w = 0. (42) 

Note that it is not the 'tHooft-Polyakov type, which is obtained through a spontaneous 
symmetry breaking by a Higgs field as in Refs. p5|- |2^JT9|1 , but the Wu-Yang type solution 
|37j. It is a kind of the dual solution of U(l) electrically charged one. The value of the 
magnetic charge is quantized as Q m = 1.0. 
On the event horizon, Eq. (|18D becomes 

A 3 7 3 (2A + l)+Aj(A — 2) 

- ^-{A 2 - 12A+ 10) = 0. (43) 
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(A 2 + 2A + 2J + (A - 2) 



We have two roots of <fi'± but find again a regular solution only for cf)' H = <p' + . The behavior 
of the field functions is similar to that of neutral and electrically charged solutions. We plot 
the M — th relation in Fig. 4. We discuss its properties in the next subsection together 
with the "colored" case. 



D. Dilatonic Colored Black Hole with Gauss-Bonnet Term 

Here we set the 'tHooft-Polyakov ansatz a = 0, i.e., purely magnetic YM field strength ex- 
ists. In the SU(2) Einstein- Yang-Mills system, the no-hair theorem for a spherical monopole 
or dyons was proved |3^j3^]. It states that there exists no static, spherically symmetric reg- 



ular non-Abelian black hole solution which has non-zero global Yang-Mills magnetic charge 
with or without an electric charge. Although it is not clear whether this kind of no-hair 
theorem holds in the present case, we set a = in this paper. 
We obtain the following relation on the horizon from Eq. (|T8|), 



C x 4>'l + C 2 <P' H + C 3 = 0, (44) 

where 

d = A 2 7 (l - w 2 H f (2AV + A 7 2 + l) - 2A 7 , (45) 

C 2 = -A 4 7 2 (l - ~ A(l - w 2 H f (2A 2 7 2 + 2Aj 2 + l) + 1, (46) 

C 3 = -^(1 - w 2 H f + Aj(l - w 2 H ) 2 (6A + 1) - 6A 7 . (47) 

This has again two roots <j)' ± , and we find a non-trivial solution only with the value <ft' + , as 
the same as before. The YM equation (|20|) on the event horizon becomes 



Then w' H is also determined by 4>h and wh- Now we have only one shooting parameter, wh, 
which will be fixed by an iterative integration with the boundary condition ([33|). Here we 
assume wh > without loss of generality, since the field equations are symmetric for the 
sign change of w. 
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With the above conditions, we solve Eqs. ([16]) ~ (|18D and (^) numerically and find 
discrete families of regular black hole solutions which are characterized by the node number 



n of the YM potential, just as for the colored black hole [Q. We show some of the solutions 
with n — 1 in Fig. 5. We shall call those solutions "colored" black holes as well. Since the 
YM field damps faster than ~ 1/r 2 , those black holes have no global color charge related 
to the gauge field, just like the colored black hole. Here we use double-quotation marks to 
distinct our new solutions from the original colored black hole. The dilaton field and metric 
functions are similar to those of other solutions discussed in the previous subsections. There 
is a region where the effective energy density becomes negative. 

The characteristic feature is that the YM potential w is almost scale invariant. This is 
recognized as follows. We shall normalize the variables by f# to see the scale invariance, 
i.e., f = f/f#, rh = m/ru, a = f#a. In the limit of r~n — > oo the field equations become 

f = 4^ B + o(4), (49) 

* , =i( 1 -tW +0 (fV (50) 



H , 



e~ 5 f 2 1 - 



2m 

f 



H 



O ( W J = 0, (51) 



2A y^ + e-^^^l + o(l)=0. (52) 



r 



r 



\' H / 

A prime with a hat denotes the derivative with respect to f. Eqs. ( f49D ~ fl5T| ) are decoupled 
from the YM field and are the same as those for the Einstein gravity with a massless scalar 
field. Because of the no-hair theorem, we find just a Schwarzschild solution with = 0. 
The YM field equation (|52]) is then the same as that in a fixed Schwarzschild background 
spacetime, i.e., 

f(f - 2M)w V + 2Mw' + w(l - w 2 ) = 0. (53) 

The non-Abelian YM field can have a nontrivial configuration although it makes no contri- 
bution to the black hole structure. Then the solution w = w*(r) = w* if /t~h) = w* (r/rn) 
is scale invariant. From our analysis this scale invariance is still found approximately even 
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for small black holes. The configuration of the YM potentials appears to be described by 
almost the same function of r/rn- Although the metric function 5 for a small black hole (e.g. 
Th = 1.4007) varies so rapidly that the above argument may seem to be invalid. This may 
be understood as follows. 5 is decoupled from other field functions even in the original basic 
equations. As a result the YM potential w is not affected by 5, (see Eq. (|5B|)). Furthermore 
if m(r) and e _7< ^ are almost constant, which is confirmed from our numerical solutions, then 
the equation for w(r) turns out to be Eq. (|53|) , resulting in a scale invariance of w(r). A 
similar behavior is found for large black holes in the EYMD system |23| . 



We show the M — rn relations of the neutral, the magnetically charged (Q m = 1-0) and 
the "colored" black holes with node number n = 1 in Fig. 6. There exists a lower bound for 
its mass in each branch. Since the solution does not exist for th — > 0, our result is consistent 
with j31|, which says that there is no particle- like solutions in the present system. 



Fig. 6(b) is a magnification around the critical point for the magnetically charged and 
"colored" black hole. Although the neutral and the electrically charged branches have a 
turning point, which is a critical point, the magnetically charged and "colored" branches do 
not. The lowest mass corresponds to the end point, where 0' ± coincide. 

Although we have not proved it analytically, we expect that a naked singularity appears 
at this end point as neutral solutions, because 8' for the black hole near critical mass appears 
to diverge on the horizon from Fig. 4(d), and in our numerical solutions. Hence the critical 
point coincides with the singular point S. 

IV. THERMODYNAMIC AL PROPERTY AND STABILITY 

In this section we investigate the thermodynamical properties of the dilatonic black 
holes with the GB term and then analyze their stability. We also discuss their fates via the 
Hawking evaporation process. 
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A. Temperature, Entropy and Stability 



The GM-GHS solutions and the non-Abelian black holes have an interesting thermody- 
namical property. That is, a discontinuity of the heat capacity of the GM-GHS solution 
appears, depending on the coupling constant of the dilaton field 7. Its critical value is just 
what we use in this paper [H. Similar properties were obtained in the EYMD system ||23|| . 
Our new black holes may possess similar interesting properties, which is the reason why we 
investigate the thermodynamical properties here. 

Although black hole thermodynamics in non-Einstein theories is not well studied, we can 
define the temperature and the entropy of the black holes in our model and show that they 
obey the first low of thermodynamics ||40|| . The Hawking temperature is given as 

T = -±- e -'» [1 - 2m' H ] , (54) 



for the metric (|T5|). The inverse temperature (3 = 1/T versus the gravitational mass is shown 
in Fig. 7. We show the branches of the neutral, the electrically charged (Q e = 1.0), the mag- 
netically charged (Q m = 1-0) and the "colored" black holes with n = 1. For comparison we 
also plot the GM-GHS solution, which have same temperature as Schwarzschild black holes 
in 7 = a/2 case. Comparing the branch of the neutral black holes with that of Schwarzschild 
black holes, we can guess that the GB term has the tendency to raise the temperature. This 
can be confirmed by seeing other branches. For example the dilatonic colored black holes in 
the EYMD system have a lower temperature than that of the Schwarzschild black holes in 



the mass range of Fig. 7 |23| . However the new solutions with "color charge" are under the 
branch of the Schwarzschild black holes. These behavior comes from the contribution of fh' 
term of Eq. ([54]) . fh' has the minimum value fa' = — 1/2 at the singular point in the neutral 
case. As for the electrically charged black hole we can see the effect of the GB term in from 
a different point of view. The thermodynamical properties of GM-GHS solutions changes 
drastically as the coupling constant 7 shifts across 7 = y/2. For 7 < \/2 the temperature is 
always lower than Schwarzschild black hole and it vanishes in the extreme limit, while for 
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7 > \/2 the temperature is always higher than Schwarzschild black hole and it diverges in the 
extreme limit. For 7 = y/2, which is also the applicable value in our case, the temperature 
coincides with the Schwarzschild case and it is finite even in the extreme limit. However the 
behavior of the temperature of the electrically charged solutions is similar to the GM-GHS 
solution with 7 > y/2. The critical value of 7 in the system including the GB term must be 
different from 7 = y/2 if such a value exists at all. 

The effect of the GB term also appears when we discuss the heat capacity. It is always 
negative for all new black holes found here in spite of the existence of the gauge field. Hence 
the GB term produces a large effect on the thermodynamical properties. 

We also calculate the entropy of our new solutions. Although the entropy is expressed 
by the quarter of the area of an event horizon in the Einstein theory, this relation cannot 



be applied in the non-Einstein theory. Here we adopt the entropy proposed by Wald fiO 
which originates from the Noether charge of the system. This entropy has several desirable 
properties. For example, it can be defined in a covariant way in all diffeomorphism-invariant 
theories, which includes the present model, and it obeys the first low of black hole thermo- 
dynamics for an arbitrary perturbation of a stationary black hole. 
The explicit form of the entropy is derived by 

S = -2tt / ^e^, (55) 

where £ is the event horizon 2-surface, e^ u denotes the volume element binormal to £ and 
the integral is taken with respect to the induced volume element on X. Ej^ upa is defined as 

= (56) 

where £ is a Lagrangian density. Substituting the present model ([[]) into (|55D and (0), we 
obtain a simple form 

A H I \ a' 

H 

where Ah is the area of the event horizon. This expression is also true for 7 = case, 
in which we find only a trivial Schwarzschild solution. Although the GB term is totally 
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S='-f[l + £ r e-«*), (57) 



divergent for 7 = and does not give any contributions to the field equations, the black hole 
entropy defined here has the additional constant 7ra'/2 via a surface term. Fig. 8(a) is the 
plots of S/ S'sch for new solutions. The entropy of neutral black hole is always larger than that 
of the Schwarzschild solution with 7 = 0, although the area of neutral black holes is smaller 
than that of Schwarzschild solutions with a same mass, from Fig. 2. Since ipu is negative 
for neutral solutions in Fig. 1(a), then the second term of Eq. 57 make the entropy much 
larger. Similarly, charged black holes have larger entropy than that of Reissner-Nordstrom 
solutions. We depict the M-S diagram in Fig. 8(c), which shows the same property as in 
the Einstein theory that the entropy becomes smaller when the system includes the gauge 
field. 

For the neutral and electrically charged case, we find the critical point C and have two 
different black hole solutions with the same mass. Since the entropy may be a good indicator 
for the stability pT7|— 19[| , we show the fine structure of M-S* diagram near the critical point 



in Fig. 8(b) for the neutral black hole. We subtract a linear function which passes through 
an appropriate point A and the end point S from the entropy S(M) in order to show the 
structure near the critical point clearly. Hence only the relative values on the vertical axis 
are significant. We find a cusp structure at the critical point C. The cusp reminds us of 
the catastrophe theory. Catastrophe theory is a mathematical tool established by Thorn 
to explain a variety of changes of state in nature, in particular a discontinuous change of 
states which occurs eventually in spite of gradual changes of parameters of a system f41fl . It 
is widely applied in various research fields and in particular we showed that this theory is 
applicable to the stability analysis of various types of non-Abelian black holes [|18|.|nj . 

As in the previous non-Abelian black hole case, we adopt the mass and the entropy 
of a black hole as a control parameter and a potential function, respectively. Then the 
existence of the cusp structure shows that the stability changes at this point. Then we can 
conclude that the upper branch, AC, is more stable than lower branch, CS, since the entropy 
is the potential function. As a result, the singular point S becomes unstable. Although the 
stability of the upper branch AC is not confirmed from Fig. 8(c) because the catastrophe 



theoretical method gives us only the relative stability, we expect that it is stable since no 
other solution branch exists. It is also understood from the following argument. We know 
the stability of the black hole solutions in both EMD and EYMD systems. As in Ref. |JT], if 
we have a numerical constant f3 in front of GB term, we may have another control parameter 
(3. If we solve a series of black hole solutions from (5 = to 1, we find that our present 
solutions on the AC branch must be stable. 

The two curves AC and CS do not intersect with a non-zero angle but become tangent 
to each other. Since a temperature is expressed by T = dM/dS from the first low of 
thermodynamics, the temperature is continuous at this critical point. This is consistent 
with Fig. 7. 

Finally, it must be noted, however, that it has not been proved whether the entropy 
defined by Eq. (|55|) satisfies the second low of the black hole thermodynamics. In particular 
in our case there is a region where the effective energy density becomes negative. 

B. Fate of Dilatonic Black Holes 

We find that the temperature of the black hole is finite for all mass range. Hence we 
expect that our black holes will not stop evaporating via Hawking radiation. Then we expect 
the following fates of our black holes: When the neutral or the electrically charged black 
hole reaches the critical mass C, then the system must shift to another phase which cannot 
be described as a static spherically symmetric regular black hole. For the magnetically 
charged and the "colored" black hole, the evaporation proceeds until the spacetime reaches 
the singular point S and a naked singularity will eventually appear on the hypersurface 
where the event horizon was located. 

However we have to study the evaporation process more carefully. Recall that the tem- 
perature of the extreme black hole with 7 > in the EMD system is infinite. Although 
the naive expectation is that this leads to an infinite emission rate, Holzhey and Wilczek 
showed that the potential, through which created particles travel away to infinity, grows 
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infinitely high in the extreme limit, and then it is expected that the emission rate could be 



suppressed to a finite value | i2| . However, it turns out that our numerical analysis shows 
this expectation is incorrect |}43| . In our case, if the potential barrier becomes infinitely large 
at the critical point C or at the singular point S, the emission rate might be suppressed to 
zero though the temperature of the black hole remains finite. Then evaporation stops and 
the black hole cannot reach the critical or the singular point. Therefore we have to calculate 
the potential barrier for some field in the background of our new solutions. 

Here we examine a neutral massless scaler field $. It obeys the Klein-Gordon equation; 



$ 'f = 0. 



The scalar field is expanded in harmonics. We study one mode of 



-Y, 



lm \ 



-iujt 



(58) 



(59) 



Eq. fl58|) becomes separable, and then the radial equation can be written as 

d 2 



dr 



*2 



+ io 2 -V 2 (r) 



X(r*) = 0, 



where r* is the tortoise coordinate defined as 



and V 2 (r) is the potential; 

V\r) = e — 1 



d 

dr* 



2m\ 
r J 



2m 



-6 



d_ 

dr 



l { l + l) + 2 -^-2m'-r5'(l-^) 
r \ r J 



(60) 



(61) 



(62) 



We calculate the potential of the new black holes. We plot the ratio of the black hole 
temperature to the maximum of the potential for the neutral solution (Fig. 9). We also plot 
the Schwarzschild black hole case for comparison. All plots for the neutral solutions have 
larger values than that of Schwarzschild black holes. This result means that the evaporation 
process of the neutral solution may be faster than in the Schwarzschild case. Of course, the 
width of the potential is another factor which determines the evaporation rate. However 
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it does not differ much from that of the Schwarzschild black hole (see Fig. 10). Hence we 
conclude that the black hole continues to lose its mass, even when the solution approaches 
the critical mass point or the singular solution one. We then may be faced with a naked 
singularity, or a transition to a new black hole state (in the neutral or the electrically charged 
case). Since we have no black hole solutions below the critical mass, the new state might 
be a time dependent evaporating black hole or just a naked singularity, which was hidden 
behind the event horizon. It is an open question. 

V. CONCLUDING REMARKS 

We have studied dilatonic black holes with a GB term and found four new types of solu- 
tions, i.e., the neutral, the electrically charged, the magnetically charged and the "colored" 
black holes. The structures of these black holes and the field configurations are almost the 
same. This may be because the GB term becomes dominant in the basic equations. As a 
result, the following previously-unknown properties are found with the GB term: 

(1) The effective energy density becomes negative near the event horizon, which may be 
responsible for the existence of a new type of black hole. 

(2) The neutral or the electrically charged black holes have the critical point C, below which 
mass no solution exists, and the singular point S, where a naked singularity appears. 

(3) For the magnetically charged or the "colored" black holes there is no critical point but 
only the singular point S where the black hole has a finite minimum mass. 

(4) The entropy is calculated. It takes a minimum value at the critical point C for the neu- 
tral and electrically charged black holes. We also find a cusp structure in the MS diagram, 
which means that the stability changes at the critical point C. Since we do not find a cusp 
structure for magnetically charged and "colored" black holes, their stabilities do not change. 

(5) The black hole temperature is always finite even at the critical or the singular point. The 
heat capacity is always negative, like the Schwarzschild black hole, even for the electrically 
or magnetically charged case. This is because the mass of the extreme black hole without 
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the GB term, which is given by GM extr eme/ V 7 ^ 7 = g 2 Qe for a fixed charge, is much smaller 
than the lower bound for the present black hole mass. Then the charge per unit mass for 
our black holes is rather small. 

(6) From the finiteness of the temperature and the study of the effective potential for a 
massless scalar field, we may conclude that the evaporation process will not be stopped even 
at the critical point or the singular point. Hence it seems inevitable that the black hole 
shifts to another state such as a dynamical evaporation phase or that a naked singularity 
will appear. 

We should mention several further points. There may be a dependence on the metric 
frame. In this paper we have worked only in the Einstein frame. However some results may 
change if we go into the string frame since the conformal transformation includes a nontrivial 
dilaton field. Therefore we studied the system again in the string frame. For example, we 
plot the M — r H relation of the neutral black holes in the string frame in Fig. 11. It is 
known that the gravitational mass and the inertial mass are not equivalent in the string 
frame. Here we use the former in Fig. 11. We can find similar structures; the existence of 
a critical point and a singular point. However, the present critical point is not the same as 
that in the Einstein frame. Since the stability should not change in a choice of the frames, 
we expect that the entropy will not take a minimum value at this critical point and no 
cusp structure will appear in the M string -S strin g diagram. From the consistency for stability 
analysis, the entropy should be minimized at the critical point in the string frame. 

Next is that the solution curve becomes vertical at the critical point in the M — r H 
diagram. Hence the rate of change of the size of a black hole becomes infinite there, i.e., if 
a black hole emits just one particle by the evaporation process (in fact it is possible even 
at the critical point because the temperature and the effective potential are finite), the size 
of the black hole will change drastically. This means that the quasi-static approximation 
is broken and the thermodynamical approach may also break down near the critical point. 
Thus we have to study this process taking the back reaction into account. 
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Finally, our results are obtained by using the model (|l|) which includes only the leading 
terms of the expansion parameter a'. This expansion may not be valid for Planck scale black 
holes. The existence and behavior of our black holes near the critical or the singular point 
will be modified. Hence there is a possibility that the critical and singular point are removed 
by taking into account the higher or all orders in a'. It might turn out that a particle-like 
solution found by Donets and Gal'tsov for (3 < 0.37 exists in a string theory (/3 = 1), which 
could be important in discussion about the final state of a black hole evaporation or the 
information loss problem. 
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Figure Captions 

FIG. 1: We plot the behavior of (a) the dilaton field (b) the mass function and (c) the 
lapse function of the neutral black holes for three different radii of event hori- 
zon; r H /^ = 1.5704 (solid line), = 1.9784 (dotted line), = 3.0556 (dashed 
line). We find the region where the effective energy density becomes negative. 
5' seems to diverge on the event horizon as the black hole gets smaller. 

FIG. 2: (a) The mass-horizon radius diagram for the neutral black holes (Q e = 0) 
and the electrically charged black holes (Q e = 0.4, 1.0). There is an end 
point, where a naked singularity appears, for each curve. We also plot the 
Schwarzschild black hole and the GM-GHS solution for comparison, (b) and 
(c) is a magnifications around the singular point S of neutral black holes and 
electrically charged black holes respectively. We find the critical point C, below 
which no solution exists. 

FIG. 3: We plot the behavior of (a) the dilaton field (b) the mass function and (c) the 
lapse function of the electrically charged black holes with Q e = 1.0 for three 
different radii of event horizon; r H /\fa' = 1.52 (solid line), = 1.70 (dotted 
line), = 2.00 (dashed line). The qualitative behaviors of these functions are 
almost same as those of the neutral case. 

FIG. 4: We plot the behavior of (a) the dilaton field (b) the mass function and (c) the 
lapse function of the magnetically charged black holes with Q m = 1.0 for three 
different radii of event horizon; r H /\fa i = 1.3864 (solid line), = 1.8929 (dotted 
line), = 3.0101 (dashed line). The qualitative behavior of these functions is 
almost same as that of the neutral and electrically charged case. 
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FIG. 5: We plot the behavior of (a) the dilaton field (b) the mass function (c) the 
lapse function and (d) the YM potential of the "colored" black holes for three 
different radii of event horizon; ru/Vo' = 1.4007 (solid line), = 1.8964 (dotted 
line), = 3.0117 (dashed line). The behavior of the YM potential w hardly 
depends on the size of the black hole. The Yang-Mills field damps faster than 
~ 1/r 2 so "colored" black holes have no global YM charge. 

FIG. 6: (a) The mass-horizon radius diagram for the neutral , the magnetically charged 
(Qm — 1-0) and "colored" black holes (n — 1). We plot the Schwarzschild black 
holes and GM-GHS solutions for comparison. There exists a lower mass for 
each branch, (b) is a magnification around the end points of magnetically 
charged and "colored" black holes. There is no critical point but only the 
singular point, where a naked singularity appears. 

FIG. 7: The mass-temperature relations for the neutral, the electrically charged (Q e = 
1.0), the magnetically charged (Q m = 1.0) and the "colored" black hole (n = 
1). The temperature remains finite in all mass range for each branch. The 
effect of the GB term also appears in the heat capacity, which is always negative 
for each branch in spite of the existence of the gauge field. 
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(a) We plot the ratio of the entropy for the neutral, the electrically charged 
(Q e = 1.0), the magnetically charged (Q m = 1.0) and the "colored" black hole 
(n = 1) to that for the Schwarzschild black hole. Here we define the entropy 
of the Schwarzschild black hole by Eq. (0), so it contains the contributions 
of a surface term. We plot the entropy of Reissner-Nordstrom black holes by 
dashed line for comparison. The neutral and the charged solutions have larger 
entropy than Schwarzschild and Reissner-Nordstom black holes respectively. 

(b) is the entropy for these four types of solutions. The entropy becomes small 
when the system includes the gauge field. This behavior is same in Einstein 
theory, (c) is the magnification around the critical point of the neutral black 
holes. We change the vertical axis S — > S* in order to show the structure 
clearly by subtracting the linear function, which passes though the points A 
and S, from the original entropy S. Hence the absolute value have no meaning. 
We can find a cusp structure, which means that the stability change occurs at 
the critical point C. 

The ratio of the temperature to the maximum of the potential, through which 
particles travel, of the neutral black holes. The ratio of new solutions is al- 
ways larger than that of the Schwarzschild black hole. Hence the evaporation 
proceeds until the black hole reaches to the critical point. 
We plot the configurations of the potential of neutral massless scalar field 
around the neutral black holes with three different radii of the event hori- 
zon; th/Vc/ = 1.5704, = 1.9784 and = 3.0556. We also plot those of the 
Schwarzschild solution case with the same size of black hole as each neutral 
solution. The configurations are similar to those of Schwarzschild solutions for 
all horizon sizes. 

The mass-horizon radius diagram for the neutral black holes in the string 
frame. There are also both critical and singular points. This critical point, 
however, does not coincide with that in Einstein frame. 
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